Quantum interference in laser-induced nonsequential double ionization in diatomic 
molecules: the role of alignment and orbital symmetry 
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We address the influence of the orbital symmetry and of the molecular alignment with respect to 
the laser-field polarization on laser-induced nonsequential double ionization of diatomic molecules, 
in the length and velocity gauges. We work within the strong-field approximation and assume that 
the second electron is dislodged by electron-impact ionization, and also consider the classical limit 
of this model. We show that the electron-momentum distributions exhibit interference maxima 
and minima due to the electron emission at spatially separated centers. The interference patterns 
survive the integration over the transverse momenta for a small range of alignment angles, and are 
sharpest for parallel-aligned molecules. Due to the contributions of transverse-momentum compo- 
nents, these patterns become less defined as the alignment angle increases, until they disappear 
for perpendicular alignment. This behavior influences the shapes and the peaks of the electron 
momentum distributions. 
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I. INTRODUCTION 

High-order harmonic generation (HHG) or above- 
threshold ionization (ATI) owe their existence to the re- 
combination, or elastic collision, respectively, of an elec- 
tron with its parent ion or molecule Such processes 
take place within a fraction of a laser cycle, which, for 
a typical Titanium-Sapphire, high-intensity laser pulse 
is T ~ 2.7/s. Therefore, HHG and ATI occur within 
hundreds of attoscconds H. As a direct consequence, 
one may use both phenomena to retrieve information 
from a molecule, in particular the configuration of ions 
with which the active electron recombines or rescatters, 
with subfemtosecond and subangstrom precision. Con- 
crete examples include the tomographic reconstruction 
of molecular orbitals [H, and the probing of structural 
changes in molecules with attosecond precision 

In particular diatomic molecules, due to their simplic- 
ity, have attracted a great deal of attention [j| H, 0, H, 
!, (13 M IH 0, Q \m M G3 El El. For such sys- 
tems, sharp interference maxima and minima in the ATI 
or HHG spectra have been identified, which could be ex- 
plained by the fact that harmonic or photoelectron emis- 
sion took place at spatially separated sites. Therefore, 
a diatomic molecule may be viewed as the microscopic 
counterpart of a double-slit experiment [T?^ | . 

Apart from the above-mentioned phenomena, nonse- 
quential double or multiple ionization, may, in principle, 
also be employed to extract information about the struc- 
ture of a molecule. This is not surprising, since laser- 
induced recollision also plays an important role in this 
case. The main difference from the previous scenarios is 
that the returning electron rescatters inelastically with 
its parent ion, or molecule, giving part of its kinetic en- 



ergy to release other electrons. In particular, for the sim- 
plest case of nonsequential double ionization (NSDI), this 
may occur through direct, electron-impact ionization, or 
by excitation-tunneling mechanisms, in which, upon rec- 
ollision, the first electron excites the second. Thereafter, 
the second electron tunnels out from an excited state. 
If more than two electrons are involved, there may ex- 
ist more complex ionization mechanisms subsequent to 
recollision. Recently, a statistical thermalization model 
has been proposed in order to account for these mecha- 
nisms and provide a simple description of nonsequential 
multiple ionization [2fJ. 

Specifically for NSDI of molecules, early measurements 
have already revealed that the total double ionization 
yields depend on the molecular species, and that there is 
an intensity region for which the predictions of sequen- 
tial models break down [2l|. Furthermore, recent NSDI 
experiments on diatomic molecules have shown that the 
shapes of the electron momentum distributions depend 
on the symmetry of the highest occupied molecular or- 
bital. This holds even if the molecular sample is ran- 
domly aligned with resp ect to the laser-field polariza- 
tion [22|]. Indeed, in [22|], very distinct electron momen- 
tum distributions have been observed for N2 and O2, 
as functions of the electron momentum components p n u 
(n = 1,2) parallel to the laser-field polarization. For 
the former species, elongated maxima along the diago- 
nal = p2\\ have been reported, while, for O2, the 
distributions exhibit a prominent maximum in the re- 
gion of vanishing parallel momenta, and are quite broad 
along the direction pxu — — P2II + const. This has been 
confirmed by theoretical computations within a classical 
framework, which reproduced some of the differences in 
the yields. 
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Subsequently, it has been found that the peak mo- 
menta and the shape of the N 2 electron-momentum dis- 
tributions changed considerably with the alignment angle 
of the molecules, with respect to the laser-field polariza- 
tion [23|. Specifically, for parallel alignment, roughly 40% 
larger peak momenta along the diagonal pi\\ = p 2 \\ have 
been observed, as compared to the perpendicular case. 
Furthermore, for perpendicular alignment, a larger num- 
ber of events in the second and fourth quadrant of the 
momentum plane (Pi|i,P2||) nas been reported. In [23| . 
these events have been attributed to excitation-tunneling 
mechanisms. 

Despite the above-mentioned investigations, NSDI in 
molecules has been considerably less studied than HHG 
or ATI, possibly due to the fact that it is far more difficult 
to measure, or to model 24{. Solely from a theoretical 
viewpoint, even for a single atom, it is a very demanding 
task to perform a fully numerical, three-dimensional com- 
putation of NSDI differential electron momentum distri- 
butions, for the parameter range of interest. Indeed, only 
very recently, this has been achieved for Helium [2a |. 
which is the simplest species for which NSDI occurs. 
An alternative to that are semi-analytical approaches, 
mostly within the framework of the strong-field approx- 
imation (SFA). Such methods are easier to implement 
and provide a more transparent physical interpretation 
of the problem. They exhibit, however, several contra- 
dictions. A concrete example is the fact that cruder types 
of electron-electron interaction, such as a contact-type in- 
teraction at the position of the ion, lead to a better agree- 
ment with the experiments as compared to more refined 
choices, such as a long-range Coulomb-type interaction 
[H, H3|- Apart from that, it is not clear whether the in- 
terference patterns due to electron emission in spatially 
separated centers survive the integration over the trans- 
verse momentum components. This is particularly im- 
portant, as in most NSDI experiments, these quantities 
are not resolved. 

These issues add up to the already vast amount of 
open questions related to HHG and ATI in molecules. In- 
deed, even for these considerably more studied phenom- 
ena, topics such as the gauge dependence of the interfer- 
ence patterns [ill, El, EH > the role of different scattering 
or recombination scenarios M, E], El, El El , and the in- 
fluence of collective effects (28| or polarization [l|| El, [3(| 
on the electronic bound states have raised considerable 
debate. Apart from that, specifically for NSDI, it is not 
even clear if the interference patterns would survive the 
integration over the transverse electron momenta. This 
is particularly important, as, in most NSDI experiments, 
the transverse momenta of the two electrons is not re- 
solved. 

In this paper, we perform a systematic analysis 
of quantum-interference effects in NSDI of diatomic 
molecules. We work within the strong-field approxima- 
tion (SFA) and assume that the second electron is dis- 
lodged by electron-impact ionization. We also employ the 
classical limit of this model. We consider frozen nuclei, 



and the linear combination of atomic orbitals (LCAO ap- 
proximation) . This very simplified model has the main 
advantage of allowing a transparent picture of the physi- 
cal mechanisms behind the interference patterns. Specif- 
ically, we investigate the influence of the orbital sym- 
metry and of the alignment angle on the NSDI electron 
momentum distributions, and whether, within our frame- 
work, the features reported in [22| and [23| are observed. 
Furthermore, we address the question of whether well- 
defined interference patterns such as those observed in 
ATI or HHG computations may also be obtained for 
NSDI, and, if so, under which conditions. 

This paper is organized as follows. In Sec. Ill Al 
we briefly recall the expression for the NSDI transition 
amplitude, which is subsequently applied to diatomic 
molecules (Sec. IIIBI) . Thereafter, we employ this ap- 
proach to compute differential electron momentum dis- 
tributions, for angle-integrated (Sec. IIII A|) . and aligned 
molecules (Sec. IIII Bp . In the latter case, we analyze 
the dependence of the interference patterns on the align- 
ment conditions. Finally, in Sec. IIVI we state the main 
conclusions of this paper. 



II. TRANSITION AMPLITUDES 



A. General expressions 

The simplest process behind NSDI corresponds to the 
scenario in which the first electron, upon return, releases 
the second by electron-impact ionization. The pertaining 
transition amplitude, within the SFA, reads 

/OO pt p 
dt dt' d 3 kV Pnk V k0 e tS{p ' M '\ (1) 
-oo J — OO J 



with the action 



2 »c 

s( Pn ,k,t,t') = -V / 

n=l Jt 



'[Pn + Afr)] 5 



dT 



[k + ACr)] 1 



dr - E 02 t - E 01 t{2) 



and the prefactors 



v k0 =< Ht')\v\4 1] > 



(3) 



and 



^P„k=<PiW,P 2 Wl^i2|k(i),4 2) >. (4) 



Eq. fl} describes the physical process in which an elec- 
tron, initially in a bound state | 0q >, is released by tun- 
neling ionization at a time t' into a Volkov state |k(f) >. 
Subsequently, this electron propagates in the continuum 
from t' to a later time t. At this time, it is driven back by 
the field and frees a second electron, which is bound at 



3 



\4>^ >, through the interaction V12. Finally, both elec- 
trons are in Volkov states. The final electron momenta 
are described by p„(n = 1, 2). In the above-stated equa- 
tions, Eo n (n = 1,2) give the ionization potentials, and 
V the binding potential of the system in question. The 
form factors ([3|) and ((4| contain all the information about 
the binding potential, and the interaction by which the 
second electron is dislodged, respectively. 

Clearly, Vko and V^k are gauge dependent. In fact, 
in the length gauge p„ (r) = p„ + A(r) and k(r) = k + 
A(t)(t = t,t'), while in the velocity gauge p„ (r) = p„ 
and k(r) = k. This is a direct consequence of the fact that 
the gauge transformation xi^v = exp[— iA(r) • r] from 
the length to the velocity gauge causes a translation p — ► 
p— A(t) in momentum space. This is particularly critical 
for spatially extended systems, such as molecules, and 
may alter the interference patterns. A similar effect has 
also been investigated for high-order harmonic generation 
[Til [29j and above-threshold ionization [H, |3(J. These 
discrepancies can be overcome by adequately dressing the 
electronic bound states, in order to make them gauge 
equivalent [l4j . In this work, we will restrict ourselves to 
the field-undressed case. We expect, however, based on 
the results of [l4j]> that the field-dressed velocity gauge 
distributions will have very similar patterns to the field- 
undressed length gauge distributions. The same will hold 
for their dressed length-gauge and the undressed velocity- 
gauge counterparts. 



B. Diatomic molecules 

We will now consider the specific case of diatomic 
molecules. For simplicity, we will assume frozen nuclei, 
the linear combination of atomic orbitals (LCAO) ap- 
proximation, and homonuclear molecules. Explicitly, the 
molecular bound-state wave function for each electron 
reads 



4 n) (r n ) = CV - R/2) + eO« + R/2) 

where n = 1, 2, e = ±1, and Cy, = 1/^/2(1 + eS(R), with 
5(R) = / U n) (r„ - R/2)l * 4"' ( r « + R/2)d 3 r. (6) 



(n), 



The positive and negative signs for e correspond to bond- 
ing and antibonding orbitals, respectively. The binding 
potential of this molecule, as seen by each electron, is 
given by 



V(r n ) = V (r n - R/2) + V (v n + R/2), 



(J) 



where Vq corresponds to the binding potential of each 
center in the molecule. 

The above-stated assumptions lead to 



(8) 



for the bonding and antibonding cases, respectively, with 
X(k(0) = / d 3 n exp[ik(t') ■ r^Voir^ (ri). (10) 



Thereby, we have neglected the integrals for which the 
binding potential Vo(r) and the bound-state wave func- 
tion 0o (r) are localized at different centers in the 
molecule. We have verified that the contributions from 
such integrals are very small for the parameter range of 
interest, as they decrease very quickly with the internu- 
clear distance. 

Eqs. (J5J) and ([9]) do not play a significant role in the 
appearance of well-defined interference patterns. This is 
due to the fact that the times t' at which the electron 
is emitted lie near the peak field of the laser field. In 
other words, the electron trajectories relevant to the mo- 
mentum distributions start near the times for which the 
electric field is maximum. For those most important tra- 
jectories, the range of k(t') is so limited that the term 
cos(k(f') • R/2) does not cross zero. In fact, we verified 
that the prefactor Vko has no influence on the interference 
patterns (not shown). 

Assuming that the electron-electron interaction de- 
pends only on the difference between the coordinates of 
both electrons, i.e., Vn = V^i^i — 1*2), one may write 
the prefactor V^k as 



- - 2C ^V 12 ( Pl - k)cos[7>(t) • R/2]^VW) 



V™ = 



P„k ( 2 „.)9/2 



(11) 



or 



r 



(a) 



p„k 



2iCV 



T V u (pi - k)sinp>(i) • R/2]^(V(t)), 



(12) 

with V(t) = Pi(t) + f>2(t) — k(t), for bonding and anti- 
bonding orbitals, respectively. Thereby, 

4 2) CP(t)) = I d 3 r 2 exp\iV(t) ■ r 2 ]^ 2) (r 2 ), (13) 



and 



V 12 (px -k) = J d 3 r^2(r)exp[z( Pl - k) • r], (14) 

with r = T\ — Y 2 - Specifically, in the velocity and length 
gauges, the argument in Eqs. (fTTjl . (fl"2"]) is given by 
V(t) = Pl + p 2 - k and P(t) = Pl + p 2 - k + A(t), 
respectively. 

The interference patterns studied in this work are 
caused by the pre- factors V^k- Explicitly, the two-center 
interference condition defined by V^k gives the extrema 



Pi(t) + pa(t)-k(t) -R = n7r. 



(15) 
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For symmetric highest occupied molecular orbitals, even 
and odd numbers in Eq. (|15[) denote maxima and min- 
ima, respectively, whereas in the antisymmetric case the 
situation is reversed (i.e., even and odd n give minima 
and maxima, respectively). The above-stated equation 
will be discussed in more detail in Sec. III.B. 

The structure of the highest occupied molecular orbital 
is embedded in Eqs. (|8"|)- (fT2"]) . The simplest way to pro- 
ceed is to consider these prefactors and the single-center 
action © . The multiple integral in dXJ) will be solved us- 
ing saddle-point methods. For that purpose, we must find 
the coordinates (t s ,t' s ,k s ) for which S(p n ,k, t, t') is sta- 
tionary, i.e., for which the conditions d t S(p n , k, t, t') = 
dfS(p n , k, t, t') — and <9k5(p„, k, t, t') = are satis- 
fied. This leads to the equations 



and 



E 



[k + A(i')] 2 = -2£oi, 



J dr[k + A(t)] = 0, 



[ P „+A(i)] 2 [k+A(*)r 



(16) 



(17) 



Em. (18) 



Eq. ([16]) gives the conservation of energy at the time 
at which the first electron reaches the continuum by 
tunneling ionization. As a consequence of the fact that 
tunneling has no classical counterpart, this equation pos- 
sesses no real solution. In the limit Eqi — > 0, the con- 
servation of energy for a classical particle reaching the 
continuum with vanishing drift velocity is obtained. Eq. 
(|17p constrains the intermediate momentum of the first 
electron, so that it returns to the site of its release, which 
lies at the geometric center of the molecule. In this spe- 
cific case, this means the origin of the coordinate system. 
Finally, Eq. (fT8| expresses the conservation of energy at 
a later time t, when the first electron rescatters inelasti- 
cally with its parent ion, giving part of its kinetic energy 
upon return to overcome the second ionization potential 
Eq2- Both electrons then leave immediately with final 
momenta p„. If we rewrite Eq. (fT5)) as 

2 2 

5>n|| + A{t)f = [k + A(t)f 2E 02 - J2 P«±> 

71=1 n— 1 

in terms of the electron momentum components paral- 
lel and perpendicular to the laser-field polarization, this 
yields the equations of a hypersphere in the momentum 
space, whose radius defines the region for which electron- 
impact ionization is classically allowed. If the transverse 
momentum components p 2 ^ are kept fixed, they mainly 
shift the second ionization potential towards higher val- 
ues, and, effectively decrease this region. For more de- 
tails, c.f. Ref. [H. 

Using the above-stated saddle-point equations, the 
transition amplitude is then computed by means of a uni- 
form saddle-point approximation (see [331 ] for details). A 



more rigorous approach would be to incorporate the pref- 
actors (ITT1) or (fi~2"|) in the action. This would lead to mod- 
ified saddle-point equations, in which the structure of the 
molecule, in particular scattering processes involving one 
or two centers, are taken into account. Recently, however, 
in the context of HHG, it has been verified that, unless 
the internuclear distances are of the order of the electron 
excursion amplitude, both procedures yield practically 
the same results [III [13]. Therefore, for simplicity, we 
will restrict our investigation to single-atom saddle-point 
equations (fll)|) - (fr8"|) , together with the two-center prefac- 
tors (TTTj) or (Pj) , 



III. ELECTRON MOMENTUM 
DISTRIBUTIONS 

In this section, we will compute electron momentum 
distributions, as functions of the momentum components 
(Pi\\iP2\\) parallel to the laser-field polarization. We ap- 
proximate the external laser field by a monochromatic 
wave, i.e., 



E(i) = So sin u>te x . 



(20) 



This is a reasonable approximation for pulses whose du- 
ration is of the order of ten cycles or longer (see, e.g. [HJ 
for a more detailed discussion). In particular, we will 
investigate how the symmetry of the molecular orbitals 
influence the electron momentum distributions. For a 
monochromatic driving field, these distributions read 

F(pi\\,p 2 \\) = JJd 2 Pl± d 2 p 2± \M R {p n ,t,t')+M L (p n ,t,t')\ 2 

(21) 

where M R (p n , t, t') is given by Eq. ([1]), and 
M L (p n ,t,t') = M R {~p n ,t± T/2,f ± T/2). The sub- 
scripts L and R denote the left and the right peaks in the 
electron momentum distributions, respectively. Thereby, 
we used the symmetry A(t) = ±A(i ± T/2), and inte- 
grated over the transverse momenta. We will also con- 
sider situations for which the transverse momenta are 
resolved. In this case, the integrals in (f2"Tj) are dropped. 

The above-stated distribution may also be mimicked 
employing a classical ensemble computation, in which a 
set of electrons are released with vanishing drift momen- 
tum and weighed with the quasi-static rate 



R(f) ~ IWr'exp \~2(2\E 01 \) 3/2 /(3\E(t')\) 



(22) 



Subsequently, these electrons propagate in the contin- 
uum following the classical equations of motion in the 
absence of the binding potential. Finally, some of them 
return and release a second set of electrons. Explicitly, 
this distribution is given by 



F c \p lh p n ) = d z p 1± d'p 2± F cl (p 1 ,p 2 ), (23) 
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with 



F Gi (pi,p 2 )= J dt Rit'Wp^lV^l 2 




E 02 \-E r {t)\{?4) 



where E T (t) = [k + A(t)] 2 /2 is the kinetic energy of 
the first electron upon return (see [26[ for details). One 
should note that the argument in Eq. (f2"l)) is just Eq. 
(I18p . which expresses conservation of energy following 
rescatter. This argument implicitly depends on t' , since 
both start and return times are inter-related. If the laser- 
field intensity is far above the threshold, i.e., if the classi- 
cally allowed region is large, both approaches yield very 
similar results 13511. 
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A. Angle-integrated distributions 

As a first step, we will discuss angle-integrated elec- 
tron momentum distributions from Eq. (fT]), for different 
gauges and orbital symmetry. To first approximation, we 
will assume that the second electron is dislodged by the 
contact-type interaction 



V 12 (r 1 - r 2 ) = 5(rj - r 2 ), 



(25) 



and that the electrons are bound in Is states. These 
assumptions have been employed in [22| , and led to a 
reasonable degree of agreement with the experimental 
data. In this case, the prefactor Vl2(Pi — k) = const, in 
(fTTj) - (|12p . and the Fourier transform of the initial wave 
function of the second electron reads 



1 



[2E 02 + V{t) 



212 ' 



(26) 



The prefactors Vko and V^> n k agree with the results in 
[22| . for which the velocity gauge was taken. 

We will consider the ionization potentials and equilib- 
rium internuclear distance of N 2l and laser-field intensi- 
ties well within the experimental range. To first approxi- 
mation, we will model the highest-occupied molecular or- 
bital of N 2 using the symmetric prefactor (|lip . In order 
to facilitate a direct comparison, we will also include the 
antisymmetric prefactor (|12p . and the single-atom case, 

( 2) 

for which V Pn k ~ Vi 2 (pi — k)<^Q ' (V(t j), and employ the 
same molecular and field parameters for all cases. 

Figure [1] depicts the above-mentioned distributions. In 
general, even though different gauges and orbital sym- 
metry lead to very distinct prefactors, the shapes of the 
distributions are very similar. This is due to the fact 
that the momentum region for which the transition am- 
plitude JT]) has a classical counterpart is relatively small. 
Indeed, we have verified that, for vanishing transverse 
momenta pi± = p 2 ± = 0, this region starts slightly be- 
low ±\/{7p, and extends to almost ±3\/ 



. This is the 

case for which the classically allowed region is the most 



FIG. 1: Angle-integrated electron momentum distributions as 
functions of the momentum components (pi||,P2||) parallel to 
the laser-field polarization, computed using the contact-type 
interaction (|25p . The field intensity and frequency have been 
taken as / = 1.5 x 10 14 W/cm 2 , and uj — 0.057 a.u., respec- 
tively, and the ionization potentials E01 = 0.573 a.u. and 
E02 = 0.997 a.u. correspond to A^2 at the equilibrium inter- 
nuclear distance R = 2.068 a.u. The upper and lower panels 
have been calculated in the velocity and the length gauge, re- 
spectively. Panels (a), and (d) correspond to the single atom 
case, panels (b) and (e) to the symmetric prefactors <(Sj and 
(|lip . and panels (c) and (f) to the antisymmetric prefactors 
© and (fit)) . 



Up the contributions to the 



extensive, so that below ± 
yield are negligible. Hence, the maxima and the shapes 
of these distributions are determined by the interplay be- 
tween phase-space effects and the prefactor (|26l) . 



In the length gauge, Eq. (|26| is very large near Pi\\ + 
p 2 |j = ±l.hyJTJp~, while in the velocity gauge this holds 
for P\\\+ P2\\ = ±0.5-v/C/p. This is in agreement with the 



±0.5 V 

features displayed in Fig. 1. In fact, a closer inspection 
of the length-gauge distributions shows that they exhibit 
slightly larger maxima, near = p 2 \\ = ±1.5-^/17^, and 
are broader along = — p 2 \\ than their velocity-gauge 
counterparts. In the velocity gauge, since the peak of the 
prefactor lies outside the classically allowed region, we ex- 
pect that the yield will be maximal near the smallest mo- 
mentum values which have a classical counterpart. This 
agrees with Figs. Ha)-(c), which exhibit peaks slightly 
above ±yJTjp. 

In Fig. [U one also notices that the distributions are 
nearly identical in the single-atom and molecular case. 
This is possibly due to the fact that the distributions are 
being angle-integrated. Apart from that, we have veri- 
fied that, within the classically allowed region, there is at 
most a single interference minimum. This may addition- 
ally contribute for the lack of well-defined interference 
patterns. 
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B. Interference effects 

For the above-stated reasons, in order to investigate 
whether interference patterns are present in the NSDI 
electron momentum distributions, we will proceed in 
many ways. First, we will increase the classically al- 
lowed momentum region, and hence the radius of the 
hypersphere given by Eq. (|19| . For that purpose, we 
will increase the intensity of the driving laser field. Sec- 
ond, in this section, we will consider aligned molecules, 
as it is not clear whether integrating over the alignment 
angle washes the interference patterns out. One should 
note that, for the parameters considered in this work, the 
De Broglie wavelength of the returning electron is much 
larger than the equilibrium intcrnuclcar distance of N 2 . 

Finally, in order to disentangle the influence of the 
prefactor which accounts for the two-center interference 

(2) 

from that of tp^ (V(t)), we make the further assumption 
that V12 is placed at the position of the ions. Without 
this assumption, prefactor (126p corresponding to the con- 
tact interaction depends on the final electron momenta, 
and thus introduces a bias in the distributions. This may 
obscure any effects caused solely by the molecular pref- 
actors. 

Explicitly, this reads 

V 12 = *(n - r 2 ) [<5(r 2 - R/2) + 5(r 2 + R/2)] . (27) 

Such an interaction has been successfully employed in 
the single-atom case, and led to "balloon-shaped" dis- 
tributions peaked near = p 2 \\ = ±2 -^/TTp. Such dis- 
tributions exhibited a reasonable degree of agreement 

with the experiments [2(|. This choice of V12 yields 

(2) 

f o = const, in addition to Vi2(pi — k) = const. 

Hence, apart from effects caused by the integration over 
momentum space, the shape of the distributions will be 
mainly determined by the cosine or sine factor in Eqs. 
(fTTj) or (|12p . The former and the latter case correspond to 
the symmetric or antisymmetric case, respectively. The 
explicit interference maxima and minima are given by 
Eq. G5)| 

We will now perform a more detailed analysis of such 
interference condition. In terms of the momentum com- 
ponents pjii , or pi± (i = 1,2), parallel or perpendicular to 
the laser-field polarization, this condition may be written 
as cos [(R/2] or sin [(R/2] , in terms of the argument £. 
Explicitly, this argument is given by 



with 



C11 



C = C11 + Ci 



cos 9 



and 



C± = pu_ sin 9 cos tp + p 2 ± sin 9 cos(tp - 



(28) 



(29) 



(30) 



In the above-stated equations, 9 gives the alignment an- 
gle of the molecule, tp corresponds to the angle between 
the perpendicular momentum pi^ and the polarization 
plane, and a yields the angle between both perpendic- 
ular momentum components. Since we are dealing with 
non-resolved transverse momenta, we integrate over the 
latter two angles. In the velocity and in the length gauge, 
k(t) = k and k(t) = k — A(t), respectively. Interference 
extrema will then be given by the condition 



(( x + ( {] )R = mr. 



(31) 



For a symmetric linear combination of atomic orbitals, 
even and odd n correspond to interference maxima and 
minima, respectively, whereas, in the antisymmetric case, 
this condition is reversed. 

An inspection of Eqs. (f2"9")l and (|3"0|) . together with the 
above-stated condition, provides an intuitive picture of 
how the interference patterns change with the alignment 
angle 9. For parallel alignment, the only contributions to 
such patterns will be due to £|| . In this particular case, 
the interference condition may be written as 



Pi\ 



■P2\\ 



Rcosi 



k(t), 



(32) 



where cos 9=1. Eq. (f3"2"|) implies the existence of well- 
defined interference maxima or minima, which, to first 
approximation, are parallel to the anti-diagonal pin = 
— p2\\. This is only an approximate picture, as k, accord- 
ing to the saddle-point equation (fT?)) , is dependent on 
the start time t' and on the return time t. Furthermore, 
since t' and t also depend on the transverse momenta 
of the electrons (see [34| for a more detailed discussion) , 
Eq. (|32[) is influenced by such momenta. Finally, in the 
length gauge, there is an additional time dependence via 
the vector potential A(t) at the instant of rescattering. 

As the alignment angle increases, the contributions 
from the term £_l related to the transverse momenta start 
to play an increasingly important role in determining the 
interference conditions. The main effect such contribu- 
tions have is to weaken the fringes defined by Eq. (|32p . 
until, for perpendicular alignment, the fringes completely 
vanish and the electron momentum distributions resem- 
ble those obtained for a single atom. This can be readily 
seen if we consider the interference condition for 9 = tt/2, 
which is 



Pl±_ COS tp + p 2 ± COS(ip + «) = —. 

H 



(33) 



Eq. (|33|) gives interference conditions which do not de- 
pend on k(t), and which vary with the angles tp and a. 
As one integrates over the latter parameters, which is the 
procedure adopted for distributions with non-resolved 
transverse momentum, any structure which may exist in 
Eq. (f3"3"| is washed out. 

In Fig. [2j we display electron momentum distributions 
computed in the velocity gauge for a symmetric high- 
est occupied molecular orbital and various alignment an- 
gles. The symmetric case is of particular interest, since, 
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FIG. 2: Electron momentum distributions as functions of the 
parallel momenta (pi||,P2||), for several alignment angles. We 
consider the velocity gauge, symmetric orbitals, and driving- 
field intensity / = 5 x 10 14 W/cm 2 . The remaining field and 
molecular parameters are the same as in the previous fig- 
ure. The position of the interference minima, estimated by 
assuming that the first electron returns at a field crossing, are 
indicated by the lines in the figure. Panel (a), (b), (c) and 
(d) correspond to alignment angles 8 — 0,8 = 30°, 8 — 60° 
and 8 = 90°, respectively. 



recently, NSDI electron momentum distributions have 
been measured for aligned N 2 molecules 23]. For par- 
allel alignment, interference fringes parallel to the anti- 
diagonal pin = —p2\\ can be clearly seen, according to 
Eq. (f3"2")) . For small alignment angles, such as that in 
Fig. Hth) , the maxima and minima start to move towards 
larger parallel momenta. Furthermore, there exists an in- 
crease in the momentum difference between consecutive 
maxima or minima, and the interference fringes become 
less defined. This is due to the fact that the term 
which washes out the interference patterns, is getting in- 
creasingly prominent. For large alignment angles, such 
as that in Fig. djc), the contributions from this term are 
very prominent and have practically washed out the two- 
center interference. Finally, for perpendicular alignment, 
the distributions resemble very much those obtained for 
the single-atom case, i.e., circular distributions peaked at 
PiU =p 2 || = ±2y% (c.f. Refs. [H|22| for details )- This 
is expected, since the term responsible for the two-center 
interference fringes is vanishing for 6 = 90°. 

The fringes in Fig. [5] exhibit a very good qualitative 



agreement with the interference conditions derived in this 
section. Furthermore, the figure shows that, for some 
alignment angles, the patterns caused by the two-center 
interference survive the integration over the transverse 
momentum components. It is not clear, however, how 
well the position of the fringes agree with Eq. (132|) quan- 
titatively, and if it is possible to provide simple estimates 
for these maxima and minima. Apart from that, it is not 
an obvious fact that the patterns survive the integration 
over the transverse momentum, and one should under- 
stand why this happens. 

In particular, the role of the intermediate momentum 
of the first electron will be analyzed subsequently. Ac- 
cording to the return condition (|17[) . this quantity de- 
pends on the start and return times of the first electron. 
Furthermore, in the length gauge, the interference con- 
dition also depends on the vector potential A{t\) at the 
return time of the first electron. For each pair (pi||,P2||), 
the emission and return times are strongly dependent on 
the transverse momenta 3J|. Apart from that, physi- 



cally, there are several orbits along which the first elec- 
tron may return, which occur in pairs. Hence, there exist 
several possible values for k. In practice, only the two 
shortest orbits contribute significantly to the yield. The 
contributions from the remaining pairs are strongly sup- 
pressed due to wave-packet spreading. However, this still 
means that the intermediate momentum, and therefore 
the position of the maxima and minima, has two possible 
values, which depend on the start and return times, and 
also on the final momentum components. 

We have made a rough estimate of the position of these 
patterns for parallel alignment, in the velocity and length 
gauges, along the diagonal p^ = p 2 \\ = p\\. This es- 
timate is given in Table 1. For symmetric highest occu- 
pied molecular orbitals, the even and odd numbers denote 
maxima and minima, respectively, while for antisymmet- 
ric orbits this role is reversed. Thereby, we assumed that 
the first electron left at peak field and returned at a field 
crossing. This gives \k\ ~ y/Up/ (0.7 '5n) in the saddle- 
point equation (|17| . Furthermore, in the length-gauge 
estimate, we took \A(t)\ ~ 2y/U^. We have verified that 
both quantities are negative for the orbits in question. 

These estimates agree reasonably well with the electron 
momentum distributions along pin — p 2 \\ = p\\. These 
distributions are depicted in Fig. [3] for several alignment 
angles, the velocity gauge, and symmetric highest occu- 
pied molecular orbitals. The positions of the minima, for 
each angle, are indicated in the figure. These minima 
have been computed employing Eq. (|32p and the above- 
stated estimate for k. For parallel alignment [Fig. [^a)], 
the position of the extrema agree relatively well with Ta- 
ble 1. This suggests that the intermediate momentum 
of the first electron, upon return, can be approximated 
by its value at the field crossing. As the alignment angle 
increases, the patterns become increasingly blurred un- 
til they are eventually washed out by the contributions 
of (±. For instance, for 9 = 30° [Fig. (2Ub)], one may 
still identify a change of slope in the distributions, at 



8 




-4 -2 2 4 



P,/[U D ] 1/2 



-4 -2 2 4 

P,/[U P ] 1 ' 2 



FIG. 3: Electron momentum distributions for parallel mo- 
menta = p 2 || = Plli non-resolved transverse momenta and 
several alignment angles. We consider the velocity gauge, 
symmetric orbitals, and the same molecule and field parame- 
ters as in the previous figure. The position of the interference 
minima, estimated by assuming that the first electron returns 
at a field crossing, are indicated by the vertical lines in the 
figure. Panel (a), (b), (c) and (d) correspond to alignment 
angles 8 = 0, 8 = 30°, 9 = 60° and 8 = 90°, respectively. For 
comparison, the yield for 8 = are indicated as the dashed 
lines in the figure. To facilitate the comparison, the yields 
have been normalized to the same peak values. 



the momentum for which the minima n = 3 is expected 
to occur. For 9 = 60°, however, the term (± has al- 
ready washed out the interference patterns. Indeed, in 
Fig. EJc), there is no evidence of interference patterns. 
Finally, for perpendicular alignment, the distributions re- 
semble very much those obtained in the single-atom case, 
as shown in [Fig. EDJd)]. 



In order to investigate the behavior of the interme- 
diate momentum k with respect to p n ±{n = 1,2), we 
will compute electron momentum distributions keeping 
the absolute values of the transverse momenta fixed. For 
simplicity, we will take 9 — and parallel momenta along 
the diagonal, i.e., pi\\ = P2\\ = p\\ • These distributions are 
displayed in Fig. 2] In this case, there exists a region of 
parallel momenta for which the yield is oscillating, be- 
tween a maximum and a minimum parallel momentum. 



Extrema 
Order n 


Parallel momentumpn/fC/p] 1 ' 2 
velocity gauge length gauge 


i 


0.513 


1.513 


2 


1.239 


2.239 


3 


1.964 


2.964 


4 


2.689 


3.689 



TABLE I: Electron momenta corresponding to the interfer- 
ence maxima and minima given by Eq. (26), in the velocity 
and length gauges, for a parallel-aligned molecule, for the 
same field and molecule parameters as in Fig. 3. The par- 
allel momenta p« have been taken to be along the diagonal 
Pill = P211 in the momentum plane, and the transverse mo- 
menta are assumed to be vanishing. If the highest occupied 
molecular orbital is approximated by a symmetric combina- 
tion of atomic orbitals, the maxima and minima are denoted 
by even and odd number, while in the antisymmetric case, this 
role is reversed, i.e., odd and even numbers denote maxima 
and minima, respectively. 



These oscillations are due to the quantum interference 
between the two shortest possible orbits along which the 
first electron may return. These orbits constitute the pair 
that has been employed in the computations performed 
in this work. The larger the transverse momenta are, the 
less extensive this region is. This is expected according 
to Eq. (fT9| ). which delimits this region (for details see 
Ref. [23 ). 

Apart from these oscillations, Fig. [4] also exhibits the 
maxima and minima caused by the spatial two-center in- 
terference. The figure shows that the position of such 
patterns is very robust with respect to the choice of p± n , 
n = 1,2. Indeed, both maxima and minima remain at 
practically the same positions, if different transverse mo- 
menta are taken. For this reason, such patterns survive 
if one integrates over the transverse momenta. In con- 
trast, the oscillations due to the temporal interference 
get washed out. For the parameters employed in the fig- 
ure, we have verified a reasonable agreement between the 
second minimum and Table 1. The first minimum is to 
a large extent washed out by the contributions of the 
events displaced by a half-cycle, i.e., which are related to 
the transition amplitude Ml,. 

Interference fringes parallel to piii = — P2II are a l so 
present in the length gauge, and for antisymmetric or- 
bitals. This is shown in the upper panels of Fig. [5J for 
parallel alignment angle. In fact, the main difference as 
compared to the symmetric, velocity-gauge case, is the 
position of such patterns, in agreement with Eq. (|32[) . 
There is also some blurring in the patterns, in the length 
gauge, possibly caused by the fact that the vector poten- 
tial A(t) depends on the return time t. This latter quan- 
tity is different for different transverse momenta. The 
patterns, however, can be also clearly identified in this 
gauge. In all cases, however, there is no evidence of a 
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straightforward connection between an enhancement or 
suppression of the yield in the low-momentum region and 
the symmetry of the orbital. For instance, in the velocity 
gauge, the yield is enhanced if the orbital is antisymmet- 
ric. The length-gauge distributions, on the other hand, 
exhibit a suppression in that region regardless of the or- 
bital symmetry. 

In the lower panels of Fig. we display the distribu- 
tions along = p 2 \\ = p\\. Similarly to the velocity- 
gauge, symmetric case, the minima and maxima of the 
distributions roughly agree with Table 1. In fact, the 
even numbers in this table roughly give the position of 
the minima in Figs. [DJe) and (f), which correspond to 
antisymmetric orbitals, while the odd numbers approx- 
imately yield the minima in Fig. Od), which display 
the length-gauge, symmetric case. Specifically for the 
length-gauge distributions [Figs. EJd) and (e)], there is 
an overall displacement of roughly 2 \fU~ v in the position 
of the patterns. This is consistent with the modified in- 
terference conditions in this case. 



C. The classical limit 

In the following, we perform a comparison between 
the S-Matrix computation and its classical limit. In the 
single-atom case, both computations led to very similar 
results, unless the driving-field intensity is close to the 
threshold intensity 35]. At this intensity, the kinetic 
energy upon return is just enough to make the second 
electron overcome the ionization potential. Therefore, 
since the intensity used in most figures is far above the 
threshold intensity, one would expect similar results. 

In Fig. [6j we display differential momentum distri- 
butions as functions of the parallel momentum compo- 
nents, computed employing the classical model. This is 
the classical counterpart of Fig. [2| in which the quantum 
mechanical distributions are depicted for the same pa- 
rameters. Indeed, for all alignment angles depicted, the 
classical and quantum-mechanical distributions look very 
similar. Hence, even though the two-center interference 
is an intrinsically quantum mechanical effect, it can be 
mimicked to a very large extent within a classical model. 
There is also a good quantitative agreement between the 
positions of the minima and maxima in both classical and 
quantum mechanical cases. This is shown in Fig. for 
parallel momenta p\\\ = P2|| = V\\i an d several alignment 
angles. For 9 = and 9 = 30° [Figs. Efa) andEJb), 
respectively], the maxima and minima agree very well 
with those in Fig. [3] The main difference, with regard 
to the quantum-mechanical case, is that, for large align- 
ment angles, the classical distributions are more localized 
than their quantum-mechanical counterparts, especially 
in the low momentum regions. For instance, in Fig. EJd), 
the yield is much lower near p|| = 0, as compared to the 
outcome of the S-Matrix computation [Fig. E^d)]. This 
discrepancy is possibly due to the fact that the classical 
model underestimates contributions to the yield near the 
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FIG. 4: Electron momentum distributions for resolved trans- 
verse momenta, as functions of the parallel momentum pjii = 
Pi\\ = P\\, for alignment angle 9 — 0. We consider the velocity 
gauge, symmetric orbitals, and the same molecule and field 
parameters as in the previous figure. For comparison, the 
corresponding single-atom distributions are presented as the 
dashed lines in the figure. The interference minima according 
to Table 1 are indicated by the vertical lines in the figure. 
The numbers in the figure indicate the transverse momentum 
components (pi±, P2±) in units of \/Vp. 



boundary of the classically allowed region. 



IV. CONCLUSIONS 

In this work, we addressed two aspects of non- 
sequential double ionization of diatomic molecules: the 
influence of the symmetry of the highest occupied molec- 
ular orbital, and of the alignment angle, on the differen- 
tial electron momentum distributions. We considered the 
physical mechanism of electron-impact ionization, within 
the strong-field approximation, and very simple models 
for the highest occupied molecular orbitals, within the 
LCAO and frozen nuclei approximations. 

For angle-integrated electron momentum distributions, 
we have shown that, for driving-field intensities within 
the tunneling regime and compatible with existing ex- 
periments [22j |. the distributions computed with sym- 
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FIG. 5: Electron momentum distributions for a parallel- 
aligned molecule (8 = 0), different orbital symmetries and 
gauges. The upper and lower panels give the contour plots 
as functions of the parallel momenta, and the distributions 
along pi|i = p 2 || = P\\, respectively. We integrate over the 
transverse momenta, and employ the same molecule and field 
parameters as in the previous figures. The interference min- 
ima according to Table 1 are indicated by the vertical lines 
in the figure. Panel (a) and (d), (b) and (e), and (c) and 
(f) correspond to symmetric orbitals in the length gauge, an- 
tisymmetric orbitals in the length gauge and antisymmetric 
orbitals in the velocity gauge, respectively. For panels (d), 
(e) and (f), the units in the vertical axis have been chosen 
so that their upper values are unity (the original values have 
been divided by 0.016, 0.01 and 0.04, respectively). 



metric and antisymmetric orbitals (prefactors (llip and 
(fl~2|) . respectively), or different gauges, look practically 
identical. This is due to the fact that, if only electron- 
impact ionization is taken into account, the momentum 
region for which this process has a classical counterpart 
is too small to allow the corresponding pre-factors to 
have a significant influence. At first sight, this is in con- 
tradiction with the experimental findings and computa- 
tions in (2^ | . Therein, a broadening parallel the anti- 
diagonal direction has been reported only for the anti- 
symmetric case, while, for a symmetric combination of 
atomic orbitals, an elongation in the direction pxil = P2\\ 
has been observed. One should note, however, that, in 
[23 | . an effective, time-dependent second ionization po- 
tential E m {t) = E 02 - 2y/2\E(t)\ is used [H]. This 
feature has not been incorporated in the present com- 
putations. It has the effect of increasing the classically 
allowed momentum region and introducing an additional 
time dependence in the prefactors and the action. 

We have also made a detailed assessment of the in- 
terference effects due to the fact that electron emission 
may occur from two spatially separated centers. In order 
to disentangle the interference effects from those caused 
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FIG. 6: Electron momentum distributions for symmetric 
highest occupied molecular orbitals and several alignment an- 
gles, as functions of the parallel momentum components (piy , 
Pa II ), computed in the velocity gauge using the classical model 
for the same field and molecular parameters as in Fig. [2] Pan- 
els (a), (b), (c) and (d) correspond to 8 = 0, 9 = 30°, 6 = 60°, 
and 9 = 90 , respectively. 



(2) 

by the prefactor ip y (V(t)), we assumed that the second 
electron was dislodged by a contact-type interaction at 
the position of the ions. We have observed interference 
fringes in the electron momentum distributions, along 
Pin = — Pall + const for all gauges and orbital symme- 
tries. These fringes are most pronounced if the molecule 
is aligned parallel to the laser-field polarization. As the 
alignment angle increases, it gets washed out by the term 
(|33[) . which, for angle- integrated momenta, is essentially 
isotropic in the perpendicular momentum plane. Con- 
sequently, the peaks of the distributions shift towards 
higher momenta, and their shapes resemble more and 
more those obtained for the same type of interaction in 
the single-atom case. We have also found that the promi- 
nence of such peaks will depend on the integration over 
the electron transverse momenta, so that some maxima 
may be more prominent than others. 

Interestingly, we are able to observe changes in the 
peak momenta of the distributions, as the alignment 
is varied, even if a single physical mechanism, namely 
electron-impact ionization, is considered. These changes 
are caused by the two-center interference effects. This 
complements recent results, in which different types of 
collisions and double-ionization mechanisms are associ- 
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ated with changes in the peaks of NSDI distributions, 
within the context of molecules [23l.[37|. Finally, for laser- 
field intensities within the tunneling regime, the distribu- 
tions obtained including only electron-impact ionization 
are far more localized than those reported experimen- 
tally, and the differences between different gauges and 
orbital symmetries are barely noticeable. In order to as- 
sess such effects, it was necessary to consider much higher 
intensities, for which other physical mechanisms, such as 
multiple electron recollisions, would also be expected to 
play a role (37[. These discrepancies may be due to the 
fact that we are not including the physical mechanism 
in which the first electron, upon return, promotes the 
second electron to an excited state, from which it subse- 
quently tunnels out. 



FIG. 7: Electron momentum distributions for symmetric 
highest occupied molecular orbitals and several alignment an- 
gles, along p!|| = p 2 \\ = p||, computed in the velocity gauge 
using the classical model for the same field and molecular pa- 
rameters as in Fig. [3] Panels (a), (b), (c) and (d) correspond 
to 9 = 0, 9 = 30°, 9 = 60°, and 9 = 90°, respectively. 
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